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We review the basic ideas on the origin of cosmic rays with energy in excess of 
^ 10^® eV, in the hght of the most recent observational findings. The limited 
statistics of events detected by the two largest experiments currently operating 
does not allow as yet to claim the detection of the GZK feature in the cosmic 
ray spectrum or the lack of it. Although extragalactic point sources seem to be 
preferred on the basis of the small scale anisotropics detected by AGASA, the 
possibility that UHECRs with energy in excess of 10^'' eV may be the result of the 
decay of supermassive particles is not ruled out by present data. It seems clear 
that if we want to have a clear idea of the origin of UHECRs, we have to auspicate 
the realization of larger, better experiments, such as the ongoing Auger project 
and the space-borne EUSO experiment. 



1. Introduction 

Ninety years after the discovery of cosmic rays, their origin is still sur- 
rounded by many questions still seeking an answer. May be the most re- 
silient of these puzzles is the origin of the highest energy end of the cosmic 
ray spectrum, that lies now near 10^° eV. The search for the sources of cos- 
mic rays at ZeV energies is made harder by two factors: first, it is difficult 
to envision astrophysical accelerators able to achieve such large energies ^, 
and second, the laws of physics predict a flux suppression at ~ 5 x 10^^ 
eV, that has not yet been clearly detected. This suppression, known as the 
GZK feature ^ , is expected to appear when cosmic ray protons start to feel 
the inelastic collisions with the cosmic microwave background (photopion 
production), a process that has a threshold at ^ 5 x 10^^ eV. While the 
first point is mainly a matter of theoretical investigation, the second issue 
is mainly an observational issue. The status of observations has been sum- 
marized in this conference ^, and we refer the reader to that contribution 
for details, while we concentrate here on the possible interpretations of the 
observations in terms of origin and propagation of UHECRs. 
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The paper is organized as follows: in §2 we address the issue of the 
statistical significance of present observations with AGASA and HiRes; in 
§3 we summarize some plausible astrophysical sources of UHECRs; in §4 
we describe the predictions and status of the art of top-down models of 
UHECR production. We conclude in §5. 



2. First shadows of the GZK feature? 

As summarized in Ref. ^, current observations by AGASA and HiRes 
show some discrepancy in both the absolute flux, and the spectral shape 
of the UHECRs. In particular the HiRes experiment seems to show some 
evidence for the presence of the expected GZK feature. The AGASA data 
show no sign of the GZK suppression, and are consistent with a power 
law extrapolation of the lower energy spectrum. Despite some claims of 
detection of the GZK feature ^, the investigation reported in ^ suggests that 
the situation may in fact be more subtle and that more data are needed 
to substantiate such claims. As illustrated in Fig. 1 (left panel), the two 
data sets have an offset that might be accomodated if a systematic error 
in the energy determination of ~ 15% in each experiment were present. In 
Ref. ^ the propagation of UHECRs has been followed through a numerical 
simulation, so to keep track of the uncertainties connected with the limited 
statistics of events of each experiment. At energies larger than 10^° eV, 
the discrepancy between the two experiments has been found to be at the 
level of 2.6(7 without any correction for systematics. If a systematic error 
in the energy determination of ~ 15% is assumed in each experiment, the 
discrepancy is reduced further. It is interesting to note that this systematic 
error would be perfectly compatible with the estimates of systematic errors 
published independently by the two collaborations Using a Montccarlo 
simulation, it is possible to evaluate the uncertainty in the predicted fluxes, 
as due to the intrinsic stochasticity of the photopion production process 
and to cosmic variance. If these uncertainties arc taken into account, the 
discrepancy between the two data sets is at the 1.8a level without correcting 
for the systematics, and at the 1.5cr level if such correction is taken into 
account (right panel in Fig. 1^) . It is possible that the first shadows of the 
GZK feature are appearing but any claim cither of detection of the GZK 
feature or lack of it is justified only at the ~ 2a level, and therefore far 
from being conclusive. 
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AGASA 7=2.6 
E>IO" :65l 
E>IO'^'' : 42.82 ± 6.45 
E>IQ-" : 2.25 ± 1.53 
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Figure 1. left) AGASA data (circles) and HiResI data (squares); right) Results of the 
simulation for a 15% systematic error in the energy determination of AGASA *. 

3. Astrophysical sources of UHECRs 

Astrophysical sources distributed homogeneously in the universe, or fohow- 
ing the observed density field ^ produce a spectrum of UHECRs with a 
GZK feature, more pronounced in general for steeper injection spectra ^. 
Therefore the firm detection of this feature in the spectrum may serve as 
the most convincing evidence for an extragalactic astrophysical origin of 
UHECRs. Here we briefly summarize the potential astrophysical sources of 
UHECRs and the acceleration processes that are supposed to be at work. 

3.1. Active and dead Galactic Nuclei 

The term Active Galactic Nucleus (AGN) is used here to identify a very 
broad class of objects, characterized by the presence of a central massive 
black hole, and fueled by an accretion flow. The presence of jets is probably 
a characteristic of many AGNs. When the fuel is almost exhausted, the 
AGN stops being active, leaving behind a dead quasar Cosmic rays can 
be accelerated in the central regions of AGNs, in the jets of some special 
class of AGNs and in the accretion disks of dead quasars. In the central 
regions of AGNs energy losses and diffusive confinement play together to 
limit the maximum energy to 10^ — 10^ GeV Hot spots in F-R II radio 
galaxies appear to be more promising sites for the generation of UHECRs, 
due to the lower photon density and lower values of the magnetic field in 
the acceleration region The hot spot is terminated by a strong shock at 
which particles may be accelerated diffusively reaching a maximum energy 
that can be as high as 10^^ eV or more. The spectrum of UHECRs at the 
Earth is expected to have a pronounced GZK cutoff, because the closest 
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of these sources is at redshift of ~ 0.3. The propagation of UHECRs from 
generic active galactic nuclei has been recently investigated by 

Recently it has been proposed that unipolar induction in dead quasars 
may energize charged particles to ultra high energies: in these objects en- 
ergy losses due to photopion production and curvature radiation are ex- 
pected to be less important than for active galaxies. Although dead quasars 
are not expected to be bright sources, it has been proposed that they 
may show gamma ray emission, due to curvature radiation of UHECRs 
during the acceleration process. Unfortunately, at present no detailed cal- 
culation of the spectrum of UHECRs generated by dead quasars exists in 
the literature, so that it is difficult to predict in a more quantitative way 
the spectrum of particles detected at the Earth. 

Other astrophysical objects may be suitable candidates as sources of 
UHECRs, although a specific model may not have been worked out in 
the details. Recently, evidence has been found for a correlation of the 
arrival directions of cosmic rays with energy above 2 x 10^^ eV with the 
spatial location of BL Lac objects and a possible correlation with EGRET 
sources Many of the BL Lacs are at large rcdshifts. It is worth stressing 
however that this correlation concerns only cosmic rays with energy below 

5 — 6 X 10^^ eV. At these energies the pathlength is large enough to 
allow for the arrival of particles from sources at redshift in excess of 0.1. 
Therefore the correlation does not imply any need for new physics, and does 
not suggest any plausible explanation for the events with energy larger than 
5 X 10^^ eV, whose sources remain unidentified. 

3.2. Neutron Stars 

Rapidly rotating strongly magnetized neutron stars are very efficient unipo- 
lar inductors, with an electromotive force (emf) that may reach 10^^ V 
although this potential is likely to be partially short-cut by the electron- 
positron pairs in the magnetosphere of the neutron star. The emf is avail- 
able in gap regions where the condition ^ • B = is violated The 
maximum achievable energy is typically around 10"'^^ eV mainly due 
to the limitation imposed by curvature radiation energy losses. In a 
phenomenological approach was adopted to rule in favor of acceleration of 
heavy nuclei to extremely high energies at the light cylinder of young neu- 
tron stars: it is well known that most energy lost by a spinning neutron star 
is not in the form of radiation but is rather converted into kinetic energy 
of a relativistic wind. For a Crab-like pulsar, the Lorentz factor of the wind 
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must be of the order of Tw ^ 10® — 10^ although it is not clear how this 
relativistic motion is achieved. If at the light cylinder a small amount of 
Iron nuclei or protons is present, the nuclei may acquire the same Lorentz 
factor of the wind, that for young neutron stars may well bo in the useful 
range F = 10^° — 10^^. In order to escape freely from the plerion these 
nuclei should be able to cross the ejecta without suffering spallation and 
photodisintegration, which results in a constraint on the magnetic field of 
the neutron star and on its rotation period In case of Iron nuclei, the 
origin may be galactic, although this possibility does not appear to be fa- 
vored by current data, due to the lack of evidence for anisotropy connected 
to the galactic disk. In other neutron star based models, where the acceler- 
ated particles are protons 24,25^ ^j^g sources are localized in other galaxies, 
and therefore the spectrum of UHECRs at the Earth has the usual GZK 
feature. 

3.3. Gamma Ray Bursts 

Cosmic rays may be accelerated diffusively at the relativistic shock front 
created by the relativistic fireball of a Gamma Ray Burst (ORB) lO'^e 
When the burst explodes in the interstellar medium, with magnetic field 
in the /itG range, the maximum energy of the accelerated particles is 
Emax ~ W^^eVBf^, where is the magnetic field in fiG. On the other 
hand, if for instance the GRB goes off in the relativistic wind of a neutron 
star 28,29^ then the magnetic field is expected to be larger and the maximum 
energy can be in the 10^° eV range Acceleration of UHECRs in GRBs 
may also occur due to other plasma physics inspired processes, such as the 
Wakefield acceleration proposed in 

If UHECRs are accelerated in GRBs, observations require the presence 
of an intergalactic magnetic field, in order to dilute in time the arrival of 
charged particles generated in the few GRB events occurring within ^ 100 
Mpc. Depending on the power spectrum associated with such field, the 
upper limits derived from Faraday rotation measurement are in the range 
10~^ — 10""'^'^ Gauss In a magnetic field of this order of magnitude, the 
average deflection angle of UHECRs is smaller than the angular resolution 
of present experiments (2-3 degrees), therefore clusters of events are ex- 
pected. If the observed multiplets are in fact the result of bursting sources, 
the higher energy particles should always reach the detector earlier than the 
lower energy ones. Although this condition is not satisfied by the AGASA 
clustered events, it was proposed that fluctuations may in fact invert 



February 2, 2008 3:54 WSPC/Trim Size: 9in x 6in for Proceedings 



Texas2002 



6 

the order of arrival of particles with energy. The probability to obtain the 
observed multiplets simply due to statistical fluctuations has however not 
been calculated yet. 

The energetic input of GRBs in the form of UHECRs was first investi- 
gated in 32,33,34 j£ Qj^gg follow approximately the star-formation history, 
then the local rate of bursts should be unable to account for the observed 
flux of cosmic rays with energy above 10^° eV. The arguments presented 
in were recently addressed in Ref. where the burst rate and energy 
depositions were taken from recent literature 35,36^ ^^^^^ ^^le spectra of UHE- 
CRs were calculated adopting a red-shift luminosity evolution of GRBs that 
follows Ref. The calculated energy injection rate found in through 
comparison with the HiRes data, for an injection spectrum E~^-'^, is a fac- 
tor ~ 6 larger than that available in gamma rays from GRBs (a factor 9 
larger if the AG AS A data are used), which may represent a serious problem 
for the GRB models requiring high radiative efflciency. 

4. The top-down alternative 

The idea that cosmic rays with extremely high energies may arise from 
the decay of supermassive unstable particles was proposed in the seminal 
paper of Ref. More recently it has become clear that these supermassive 
particles may be either the result of the decay or annihilation of topological 
defects or quasi-stable relics of the inflationary era 41,44,45,46,47^ rpj^^ basic 
idea, common to all top-down models, is that the decay of a supermassive 
particle results in the production of a quark-antiquark pair that hadronizes 
into mesons and protons. At the source, the composition of the produced 
particles is dominated by gamma rays and neutrinos, while only about 
5% of the energy goes into protons. After propagation over cosmological 
distances, the relative abundance of protons and gamma rays is such that 
protons dominate up to energies in excess of 10^° eV, while at even higher 
energies, the composition is expected to be gamma ray dominated (see Fig. 
2). This is the case for most topological defects models. In these models, the 
gamma ray flux is in fact uncertain mainly as a consequence of the difficulty 
in measuring the low frequency radio background from within our galaxy 
(we are screened by free-free absorption). In Fig. 2 the two sets of curves 
labelled as -y-high and j-low should bracket the range of uncertainty in the 
radio background at the relevant frequencies. Much discussion exists on 
which topological defects may generate the observed fluxes, as summarized 

in 40,38 
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log[E(eV)] E [eVJ 

Figure 2. Left) Spectra of protons and gamma rays from necklaces right) Spectra 
of gamma rays from decay of super massive relic particles in the halo 

Supermassive particles can also be produced in the early universe inde- 
pendently of topological defects, either through gravitational production 
or through direct coupling to the inflaton field. Relic particles produced 
gravitationally have mass mx < H{t) < m^, where H{t) is the Hubble con- 
stant and 777.0 is the inflaton mass; if their decay time is sufficiently long, 
they can be natural candidates for cold dark matter ^^'^^ -jj^g jQj^g (jgcay 
times may be the result of some weakly broken discrete symmetry (such 
as R-parity for the case of neutralinos). Supermassive relics accumulate in 
dark matter halos, and in particular in the halo of the Milky Way. Their 
rare decays may be sufficient to account for the observed fluxes of UHECRs 
above 10^° cV '^^ In this case the composition is expected to be dominated 
by gamma rays: this is not in contradiction with the findings of that are 
stringent for lower energies (the constraint may be important for low values 
of the mass, as for instance for the case plotted as a solid line in the right 
panel of Fig. 2). The strongest signature of the model is the anisotropy 
due to the asymmetric position of the sun in the Galaxy 50,40,5i,52 dj-rent 
observations are not yet stringent enough to rule against or in favor of the 
model. 

5. Discussion and conclusions 

The spectrum of UHECRs and in particular the presence of a GZK feature 
will be properly measured by the next generation cosmic ray experiments, 
namely the Auger project ^"^^ and the space-borne EUSO observatory 
The increase in the statistics that will be achieved by these two observa- 
tional enterprises is well illustrated in Fig. 3, where the simulations in ^ 
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have been used in order to predict the error bars expected after three years 
of operation of Auger and EUSO. The question about the origin of UHE- 
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Figure 3. Expected simulated performances of Auger (left panel) and EUSO (right 
panel)'*. 



CRs will be answered when the spectrum in the energy region above few 
10^® eV will be reliably measured: if the GZK feature is there, then we will 
have a proof that UHECRs are extragalactic protons generated by astro- 
physical sources. The high statistics of large multiplicity clusters of events 
will on the other hand give a direct hint on the type of sources that we 
have to look for, because the combination of the spectrum and small scale 
anisotropics provides a potentially powerful tool to determine not only the 
energy injected per unit volume but also the energy per source. On the 
other hand, if no evidence for the GZK feature will be found, then it will 
probably be hard to avoid to invoke some kind on new physics. In this case 
the study of composition might be the tool to discriminate among different 
models. 
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